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Using a recently developed algorithm for generic rigidity of 
two-dimensional graphs, we analyze rigidity and connectivity 
percolation transitions in two dimensions on lattices of linear 
size up to L = 4096. We compare three different universal- 
ity classes: The generic rigidity class; the connectivity class 
and; the generic "braced square net"(GBSN). We analyze the 
spanning cluster density Poo , the backbone density Pb and the 
density of dangling ends Pd ■ In the generic rigidity and con- 
nectivity cases, the load-carrying component of the spanning 
cluster, the backbone, is fractal at pc, so that the backbone 
density behaves as B ~ {p ~ Pc) for p > pc. We estimate 
l3'gr = 0.25 ± 0.02 for generic rigidity and /3^ = 0.467 ± 0.007 
for the connectivity case. We find the correlation length ex- 
ponents, Vgr = 1.16±0.03 for generic rigidity compared to the 
exact value for connectivity Uc = 4/3. In contrast the GBSN 
undergoes a first-order rigidity transition, with the backbone 
density being extensive at pc, and undergoing a jump discon- 
tinuity on reducing p across the transition. We define a model 
which tunes continuously between the GBSN and GR classes 
and show that the GR class is typical. 



I. INTRODUCTION 

Scalar Percolation is a simple model for disordered sys- 
tems, and has received much attention in the last two 
decades pj-pf. This model describes the transmission of 
a scalar conserved quantity (for example electric charge, 
or fluid mass) across a randomly diluted system. How- 
ever in the calculation of mechanical properties force ( 
i.e. a vector) must be transmitted across the system [Q. 
It was originally suggested [^ that the critical behavior 
of the elastic moduli of a percolating system should be 
equivalent to that of its conductivity, but this is only 
valid for the scalar limit of the Born model of elastic- 
ity pj, a model which is not rotationally invariant and 
in many cases inappropriate. Elastic percolation is not 
in general equivalent to scalar percolation. This was first 
made clear by the work of Feng and Sen H] , which showed 
that central-force elasticity percolation was in a different 
universality class than scalar percolation, and provided 
the starting point for a renewed interest in this problem. 

It soon became clear that the elasticity problem can 
be divided in two categories ||l^ , according to the kind of 
forces which hold the lattice together. If angular forces 
are important ||q- |l^ , [l5|| , a singly-connected path across 
the lattice is enough to ensure rigidity, so any config- 
uration of bonds which is connected is also rigid. In 
this case, the geometry of the elastic backbone is ex- 
actly the same as that of the scalar percolation prob- 
lem |,|ll],p|jl|L This is the case for bond-bending [|J|] 
and beam |10|,[15| models. Thus, the elasticity percola- 



tion problem with angular forces is now well understood 
and that understanding has borrowed much from the ge- 
ometrically equivalent scalar percolation problem. It is 
the purpose of this paper to analyze the more challeng- 
ing central-force rigidity percolation transition, stressing 
the similarities and differences between this problem and 
scalar percolation. 




FIG. 1. The six bodies shown in this figure are rigidly 
connected, i.e. they belong to the same rigid cluster. But the 
removal of any bond (thin black lines) leads to the collapse of 
the structure, which is then reduced to a collection of six in- 
dependent rigid clusters (no two are rigidly connected). This 
means that the existence of a rigid connection between for ex- 
ample clusters e and / cannot be decided on local information 
only, since it depends on the presence of 'far away' bonds, i.e. 
bonds not connected to clusters e or /. 



If rigidity I3jlj,|lj,g0|-g3|j25|j2j,|3g,|35[ is provided by 
central forces alone (e.g. rotatable springs), single- 
connectedness is not enough to ensure rigidity. In this 
case a lattice that is conducting usually would not sup- 
port an applied stress (i.e. it would not be rigid). This 
was first shown by Feng and Sen H who found that the 
rigidity threshold is significantly larger than the conduc- 
tivity threshold. An exception worth mentioning is the 
case of elastic lattices under tension, or equivalently, sys- 
tems in which all springs have zero repose length p6[ . 
For such systems, conductivity and Young modulus are 
equivalent, i.e. rigidity appears at the scalar percolation 
point. It has been recently emphasized p7| that cntropic 
effects can give rise to similar effects in central-force sys- 
tems with nonzero repose length and finite temperature, 
although the connection with ref. ^^ was not established. 

The main difficulties associated with central-force elas- 
ticity are as follows: Whereas in the scalar connectivity 
case it is a trivial problem to determine when two sites be- 
long to the same connected cluster, in the case of central- 
force rigidity it is not in general easy to decide whether 
two objects are rigidly connected. For example it takes 
some thinking to see that the six bodies in Fig. |l| form 
a rigid unit. Thus it is not easy to see how a computer 
algorithm can be devised to identify rigid clusters. 

In the scalar connectivity problem, the removal of a 



singly-connected bond leads to the separation of a con- 
nected cluster into two clusters. In the rigidity case, the 
removal of an analogous "cutting bond" may produce the 
collapse of a rigid cluster to a collection of an arbitrarily 
large number of smaller ones (we call this the house-of- 
cards mechanism). Fig. |l| shows a simple example of this 
situation. Due to this fact, the transmission of rigidity 
can be "non- local" ||l4| , |32|| , since a bond added between 
two clusters on one side of the sample may induce rigidity 
between two clusters on the other side of a sample. 





FIG. 2. Three bars are in general enough to form a rigid 
connection between two rigid bodies (case a), but for partic- 
ular, degenerate cases (case b), rigidity fails even when the 
system has the right number of bars. Case b fails to be rigid 
because the three bars happen to have a common point. A 
structure formed by two bodies connected by three bars is 
generically rigid in two dimensions if it is rigid "for most ge- 
ometrical arrangements", i.e. leaving aside degenerate config- 
urations such as b), which occur with probability zero in the 
configuration space. 

A second source of difficulties in the problem of central- 
force rigidity is the existence of particular geometrical 
arrangements for which a system may fail to be rigid 139] 
even if it is rigid for most other cases ||2^,|2^,^,Q . Take 
for example two rigid bodies connected by three bars in 
two dimensions, as shown in Fig. g. This structure is 
in general rigid, but if the three bars happen to have 
a common point, then the structure is not rigid, since 
this common point is the center of relative rotations 133] 
between the two bodies. 

Particular geometrical arrangements (such as Fig. i)), 
which are non-rigid even when the structure is rigid for 
most other configurations, are called degenerate config- 
urations. Degenerate configurations appear with proba- 
bility zero if the lattice locations are "randomly chosen" . 
A lattice is thus said to be generically rigid, if it is rigid 
for most geometrical arrangements of its sites. Generic 
rigidity only depends on the topology of the underlying 
graph, i.e. ignores the possibility of degeneracies. Since 
degenerate configurations are always possible on perfectly 
regular lattices, we will assume our lattice sites to have 
small random displacements, in which case generic rigid- 
ity applies. 

Up to recently, no simple algorithms existed for the 
determination of the rigid-cluster structure of arbitrary 
lattices. Due to this, direct solving of the elastic equa- 
tions was one of the few methods H] available to obtain 
information on the structure of rigidly connected clus- 
ters. But this is very time-consuming and did not allow 
the study of large lattices. Previous simulations were 
for example not sufficiently precise to confirm or reject 
the proposal ]I^-|20|] that bond-bending and central-force 



elastic percolation might after all still be in the same uni- 
versality class. This suggestion was not inconsistent with 
some numerical results obtained on small-sized systems 

11-11 ■ 

Recently there has been a breakthrough in the sys- 
tem sizes accessible to numerical analysis p^-|25|, fol- 
lowing the development of efficient graph-theoretic meth- 
ods for the problem or generic rigidity in two dimensions 
p^,pG,Bl|. Using such methods, we study the central- 
force rigidity percolation problem on randomly diluted 
triangular lattices of linear size up to L = 3200, and con- 
nectivity percolation and body-joint rigidity percolation 
on site-diluted square lattices of size up to L = 4096. Our 
numerical algorithm ]30| ] is complementary to the "pebble 
game" |24| , |3l| , which is an implementation of Hendrick- 
son's matching algorithm in the original "joint-bar" rep- 
resentation of the network |2^(see below). This paper is 
an elaboration and extension of our two recent letters on 
this subject p3| , p5| . We extend and elaborate upon the 
numerical data presented there in several ways: by com- 
paring rigidity and connectivity percolation, by studying 
significantly larger lattices for rigidity percolation, by giv- 
ing data on site and bond diluted lattices with a variety 
of boundary conditions; by presenting results on a new 
body-joint model which is in the universality class of bar- 
and-joint rigidity percolation and; by presenting detailed 
results for a model which continuously tunes between the 
braced square lattice (which has a first order rigidity tran- 
sition) and the isotropic triangular lattice (which has a 
second order rigidity transition). 

The numerical method is briefly described in Section O. 
In Section [II , our results are presented and their implica- 
tions discussed. A comparison is made with other avail- 
able numerical and analytical results for the central-force 
rigidity percolation problem. We also discuss the issue of 
first-order rigidity, which has been the subject of a com- 
ment and reply in physical review letters p6| . Section IV 
contains our conclusions. 



II. THE NUMERICAL METHOD 

We take an initially depleted triangular lattice and add 
bonds (in the bond-diluted case) or sites (in the site- 
diluted case) to it one at a time, and use a graph-theoretic 
matching algorithm pO| in order to identify the rigid clus- 
ters that are formed in the system. For the case of bond 
dilution, p is the density of present bonds, while in site 
dilution it indicates the density of present sites. In the 
site-diluted problem, a bond is present if the two sites it 
connects are. 

We use the body-bar version Ig^l of a recently proposed 
rigidity algorithm [ p7[ . This algorithm, being combina- 
torial in nature, allows us to identify sets of sites which 
are rigidly connected, without providing any information 
on the actual values of the stresses when an external load 
is applied. The body-bar algorithm sees the lattice as 
a collection of rigid clusters (or "bodies" ) connected by 
bars, instead of points connected by bars as proposed 
in the original algorithm [ETl and as implemented in the 
"pebble game" [gl|. The body-bar representation allows 
a more efficient use of CPU and memory, as each rigid 
cluster is represented as one object. The matching iden- 
tifies rigid clusters and condenses them to one node as 
new bonds are added to the network. 

In two-dimensional rigidity, a rigid cluster has 3 degrees 



of freedom, while a point-like joint has 2. Therefore the 
minimum number of bonds needed to rigidize n joints in 
2d is (2n — 3). Matching algorithms p7pC|j3l| ] are based 
on this sort of constraint-counting. 
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FIG. 3. Infinite percolation clusters which lie in different 
universality classes: a) Connectivity percolation {g — G = 1) 
on a triangular lattice; b) Rigidity percolation {g = 2,G — 3) 
on a triangular lattice; c) Rigidity percolation {g = 2,G — 3) 
on a braced square lattice. For a) and b), boundary conditions 
are periodic in the horizontal direction while for c) they are 
free. The system size L = 64 and rigid bus-bars are set on 
the upper and lower ends of the sample. The backbone, is 
composed of 'blobs' of internally stressed bonds (thick black 
lines), rigidly interconnected by cutting bonds (gray lines). 
Cutting bonds are also called red bonds. Removing one of 
them produces the collapse of the system. Dangling ends (thin 
lines) are rigidly connected to the backbone, but do not add 
to the ability of these networks to carry a DC external load 
(or current). 

The body-bar algorithm [Bol, can be extended to han- 



dle "rigidity problems" with arbitrary values of g (num- 
ber of degrees of freedom of a joint) and G (degrees of 
freedom of a rigid cluster). Connectivity for example, 
is just a special (simple) case of rigidity with g = I and 
G = 1: the minimum number of bonds needed to connect 
n points is n — I in any dimension. Connectivity perco- 
lation can thus be studied using this algorithm. More 
details on the application of matching algorithms for the 
specific case of connectivity percolation can be found in 
Ref. 11. 

There are several ways to define the onset of global 
rigidity in a network [^. We have used two distinct 
methods. First we determine whether an externally ap- 
plied stress can be supported by the network, which we 
call applied stress (AS) percolation. Secondly we studied 
the percolation of internally-stressed (IS) regions. 

At the AS percolation point, an applied stress is first 
able to be transmitted between the lower and upper sides 
of the sample. As we add bonds one at a time, we are 
able to exactly detect this percolation point by perform- 
ing a simple test |30| which consists in connecting an 
additional fictitious spring between the upper and lower 
sides of the system. This auxiliary spring mimics the ef- 
fect of an external load, and therefore the first time that 
a macroscopic rigid connection exists, a globally stressed 
region (the backbone) appears. 

The IS critical point is defined as the bond- or site 
density at which internal stresses percolate through the 
system. This means that the upper and lower sides of 
the system belong to the same self-stressed cluster p3| , 
and this is trivially detected within the matching algo- 
rithm [Q . The AS and IS definitions of percolation are 
in principle different, but we found p3| that the average 
percolation threshold and the critical indices coincide for 
large lattices. Similar definitions apply to the connectiv- 
ity case, with the AS case being the usual definition, i.e. 
the onset of electric conductivity, and IS being percola- 
tion of "Eddy currents" . 

We define the spanning cluster (Fig. |3|) as the set of 
bonds that are rigidly connected to both sides of the 
sample. However only a subset of these bonds carry the 
applied load. This subset is called the backbone. The 
backbone will in general include some cutting bonds, so 
named because the removal of any one of them leads to 
the loss of global load carrying capability. Cutting bonds 
attain their maximum number exactly at Pc [ p8[ . The 
backbone bonds which are not cutting bonds are parts of 
internally overconstrained blobs. In the rigidity case, the 
smallest overconstrained cluster on a triangular lattice is 
the complete hexagonal wheel (twelve bonds), while in 
the connectivity case it is a triangle (i.e. the smallest 
possible loop). The spanning cluster also contains bonds 
which are rigidly connected to both ends of the sample 
but which do not carry any of the applied load. These 
are called dangling ends. This classification is standard 
in connectivity (scalar) percolation [p7[ . 

In this work we analyze several other boundary condi- 
tions, particularly in the generic rigidity case on triangu- 
lar lattice. In that case, for site dilution we analyze: AS 
with rigid bus-bars at the ends of the sample, AS with- 
out bus bars (any-pair rigidity), and IS with bus-bars. 
For bond dilution, only the AS with bus-bars case was 
studied. We determine the exact percolation point (AS 
or IS) for each sample, so we can identify and measure 
the different components of the spanning cluster exactly 
at Pc for each sample. This should be contrasted with 



usual numerical approaches, in which averages are done 
at fixed values of p, and < Pc > is obtained from finite- 
size scaling (e.g. data collapse). In that case, it is known 
that slight differences in the estimation oi < pc > can 
lead to important deviations in critical indices ||l^. This 
source of error is absent in our measurements. Sample 
averages are done over approximately 10*/i^ samples. 



III. RESULTS 

We first analyze the size dependence of the three key 
probabilities Pg, the backbone density, Po, the dangling 
end density and Poo, the infinite cluster density, exactly 
at their percolation thresholds as described in the previous 
section. In Fig. ^-c, this data is presented for the three 
different universality classes in Figs. Sa-c. 



One clear feature of Fig. | is that the BSN (Fig. |c) has 
a qualitatively different behavior than the other cases. 
For the BSN, Pg, Poo and Pd all have a finite density 
at large L, indicating that the rigidity transition is first 
order in this case. In contrast, in both the connectiv- 
ity(Fig. |ja) and rigidity cases (Figs. 0b, d), Pb and Poo 
are decreasing in a power law fashion over the available 
size ranges. However the behavior of Pd is more complex. 
First we discuss the behavior of Pb ■ 

At a second order phase transition, finite-size scal- 
ing theory predicts Pb{pc) ~ L^^ ^^ ■ Taking into ac- 
count correction-to-scaling terms, which we assume to be 
power-law, we may generally write 



Pb = CiL-^{1 + C2L- 



(1) 



This expression is fitted to our numerical data by 
choosing the set of parameters {Ci,C2,e,U!} that min- 
imize the error 
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FIG. 4. Density of backbone bonds (circles), dan- 

gling bonds (squares) and infinite cluster bonds (diamonds) 
at the AS critical point for: a) Connectivity percolation 
(gr = G = 1) on a site-diluted square lattice; b) Rigidity 
percolation (g = 2, G = 3) on a site-diluted triangular lattice; 
c) Rigidity on a randomly braced square lattice; d) Body 
rigidity (G = g = 3) on a site-diluted square lattice . 

Case c corresponds to the generic braced square net 
(GBSN), which is a square lattice to which diagonals are 
added at random with probability p. The non-generic 
version of this problem has been studied by many au- 
thors |p3|-p5[ , and it is well known that the number of di- 
agonals needed to rigidize it is not extensive: pc '^ when 
L — > oo. This is confirmed by our numerical simulations, 
which correspond to the bus-bar boundary condition. 

In Fig. 0d we also present data for the rigidity case 
5 = 3, G = 3 on a square lattice, to further test whether 
the rigidity class is universal in two dimensions. In this 
model each site of a square lattice is a body and so has 
g = i degrees of freedom. Each of these bodies is con- 
nected to each adjacent body by two bonds or bars, i.e. 
two contiguous bodies are pinned at a common point. 
Maxwell counting Eq] then implies f — 3 — 4p, so that 
the Maxwell estimate of the bond percolation threshold 
is 3/4. Our numerical estimate is Pc = 0.74877 ±0.00005, 
thus confirming the accuracy of the Maxwell approxima- 
tion. 



A plot of —log{PB)/log{L) vs. l/log{L) should then 
have an asymptotic {L — * oo) intercept equal to the lead- 
ing exponent e. Similar fitting procedures were used to 
produce Figures and 0, where the leading exponent is 
j3 /v and 1/v respectively. 

A fit of the data in Fig. Hb,d produces a rather universal 
estimate P'gr/v = 0.22 ± 0.02. In consequence the rigid 
backbone is fractal at Pc, with a fractal dimension Db — 
1.78±0.02. In the connectivity case (Fig. |a), we find [|3| 
[3'/v = 0.350 ± 0.005, or Db = 1.650 ± 0.005, which is 
consistent with the most precise prior work |[44| , [45|] . 

Now we consider Poo and Pd ■ In the connectivity case 
(Fig. |ja), an analysis of the dangling ends and infinite 
cluster probabilities (Fig. 0a) both lead to the estimate 
P/v = 0.10 — 0.11, in agreement with the exact result 
5/48. In the rigidity case however, there are strong finite 
size effects and even at sizes of L = 3200 (joint-bar rigid- 
ity. Fig ||b), and L = 4096 (body-joint rigidity. Fig. ^, 
it looks as though the dangling probability may be satu- 
rating, while the infinite cluster density continues to de- 
crease. Since Poo = Pb + Pd, it is expected that asymp- 
totically Poo and Pd must behave in the same manner. 

Clearly the numerical results for the range of system 
sizes currently available are still controlled by finite size 
effects, and the results depend on the analysis method 
chosen. Jacobs and Thorpe ||2J| chose to interpret the 
infinite cluster probability as being key. A fit to the Poo 
data of Fig. |b,d yields (3/v = 0.147±0.005 (See Fig. |b,c) 
in agreement with Jacobs and Thorpe. But a similar fit 
of the dangling end density gives (3/v ~ 0.03 for the joint- 
bar rigidity case (Fig. ^) and (3 ^ 0.01 for the body-joint 
rigidity case (Fig. Be). In our previous work Pq | we were 
guided by the Cayley tree results [Q which indicated a 
first order jump in the infinite cluster probability. We 
thus chose to interpret Fig. |4p,d as indicating a satu- 
ration of the infinite cluster probability at the dangling 
end value of about 0.1. Having extended our data from 
L = 1024 to L ~ 4096, it now looks more likely that a 
small value of (3 occurs in the rigidity case (Fig. ^,c), 
though much larger simulation sizes are required to find 
j3 /v precisely. 

Due to the slow finite size effects found in the analysis 
of the infinite cluster and dangling end probabilities, it 
is natural to be concerned about the effect of boundary 



conditions and other, usually non-universal, parameters 
on the observed results. For generic joint-bar rigidity case 
on triangular lattices we thus tested a variety of differ- 
ent boundary conditions for both site and bond dilution. 
This data is presented in Fig. g, from which it is seen 
that the conclusions drawn from the case of rigidity per- 
colation with applied bus bars are quite robust. 
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FIG. 5. The spanning cluster density exponent jS/v as 
numerically estimated for a) connectivity percolation on a 
square lattice (7^""*^ = 4096), b) rigidity percolation on a 
triangular lattice (L™'^== — 3200) and c) body rigidity on a 
square lattice {^1/^°-^ — 4096). Two estimates result in each 
case from fitting the scaling of spanning cluster density (tri- 
angles) and dangling end density (circles). Solid lines are fits 
using Eq. (§) 

Finally, the behavior of the dangling end density as a 
fmiction of p is also quite striking. This data is presented 
in Fig. 0. At very high p, nearly all bonds belong to the 
backbone, so the dangling end density approaches zero. 
Below pc, there is no infinite cluster, so there are again 



no dangling ends. There then must be a maximum in 
the density of dangling ends between pc and p = 1. As 
seen in Fig. M, the interesting feature is the abrupt drop 
in the dangling density ai pc, a feature that appears to 
become more pronounced with increasing sample size. It 
is tempting to interpret this as definitive evidence of a 
first order rigidity transition, but it is also consistent with 
the strong finite size effects seen in Figs. ||b,d and Fig. 6, 
so we must await large lattice simulations for a definitive 
analysis. 
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FIG. 6. Backbone density Pb (solid lines) and dan- 
gling-end density Pd (dashed lines) as a function of sam- 
ple size at the percolation threshold of each sample, on tri- 
angular lattices for: bond-diluted AS with bus-bars (cir- 
cles), site-diluted AS with bus-bars (diamonds), site-diluted 
AS without bus-bars (triangles) and site-diluted IS without 
bus-bars. The AS percolation point without bus-bars is de- 
fined as the concentration of sites or bonds for which there is 
for the first time a rigid connection between at least one pair 
of points on opposite sides of the sample. 
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FIG. 7. Fraction of dangling ends on the {g — 2,G — 3) 
generic rigidity infinite cluster, as a function of p, for 
site-diluted triangular lattices of size L = 32, 64, 128, 256, 
512 and 1024. Data shown here are for the AS case with 
bus-bars. 

Now we turn to the calculation of the correlation 
length exponent for rigidity percolation. When there 
is a second-order rigidity transition, there is a di- 



verging correlation length ^ ^ |p ~ Pel 



We can 



find the exponent ly of this divergence by measuring 



the sample-to-sample fluctuations in pc as a function 
of L. The dispersion cr(L) = y^ <pI>l - <Pc >i), 
and according to finite-size-scahng a{L) ~ L^^/^ . 
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FIG. 8. The thermal exponent 1/u as numerically esti- 
mated for a) connectivity percolation [g = G = 1) on square 
lattice, b) rigidity percolation (g = 2, G = 3) on a triangular 
lattice and c) body rigidity on a square lattice (G = g = 3). 
Two independent estimates result in each case from fitting the 
scaling of red bonds (triangles) and fluctuations in pc (circles) . 

An asymptotic analysis for a{L) is shovifn in Fig. pk for 
connectivity percolation, in pp for joint-bar rigidity and 
in |8|c for body-joint rigidity. From these figures we esti- 
mate l/v — 0.75 ±0.01 (the exact value is l/v = 3/4) for 
connectivity percolation and \/v = 0.85±0.02 for rigidity 
percolation. This provides further strong evidence that 
rigidity percolation is second order in two dimensions, 
though not in the same universality class as scalar perco- 
lation. In the case of the first order rigidity on the braced 
square net (Fig. ||c), the variations in pc behave as L"'^/^, 
in accordance with analytical results for this model Go| . 
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FIG. 9. The volume fraction of a) backbone bonds b) 
cutting bonds, and c) the fluctuation o{pc) at the rigidity 
threshold for: Square lattice with g = - this is the braced 
square net (filled squares); Square lattice with g = 0.1 (cir- 
cles); Square lattice with q = 0.40 (squares); Triangular lattice 
with bond dilution (diamonds) and Triangular lattice with site 
dilution (triangles). 

Our algorithm also identifies the cutting (also called 
red or critical) bonds at the percolation point, for the 
case of AS percolation. The number Nn of red bonds 
scales at pc as L^ . Coniglio [pSl has shown that x = 
l/v exactly, for scalar percolation. Numerical evidence 
suggesting that x — \/v also in rigidity percolation was 
first presented in [g3|. It is in fact possible to extend 
Coniglio's reasoning to the case of central-force rigidity 
percolation |Q. It turns out that x = l/v has to be 
rigorously satisfied also in this case, and therefore a{L) 
and 1/Nr{L) must have the same slope in a log-log plot. 
Analysis of the number of cutting bonds is also presented 
in Fig. pi and yields values of l/v consistent with the 
analysis of variations in percolation thresholds described 




in the previous paragraph. 

Since the Cayley tree model |3q] g ives behavior quite 
similar to the braced square net ]35|], i.e. a first-order 
rigidity transition, it is interesting to ask whether the 
rigidity transition is "usually" like that on the braced 
square net (i.e. first order), or whether the second order 
transition found on triangular lattices is more typical. In 
order to probe this issue, we analyze a model which inter- 
polates between the braced square net and the triangular 
lattice. In the braced square net, the random diagonals 
are present with probability pd, to make the lattice rigid 
it is sufficient (though not necessary) to add one diagonal 
to every row of the square lattice. The probability that 
a spanning cluster exists is then P+ = (1 — (1 — pd)^)^ , 
from which we find pd* ~ InL/L. 

We generalize this model by randomly adding the di- 
agonals (with probability pd) to a square lattice whose 
bonds have been diluted with probability q. The braced 
square net is q = 0, while ii q ~ 1 — pd this model is 
equivalent to the bond-diluted triangular lattice. Typi- 
cal results for various values of q are presented in Fig. |^. 
It is seen that even for a small amount of dilution of the 
square lattice, e.g. q — 0.10, the rigidity transition re- 
turns to the behavior characteristic of the homogeneously 
diluted triangular case (see Fig. ^). We find that for suf- 
ficiently large lattice sizes, the universal behavior found 
in the other rigidity cases holds for any finite q < 0.5 (for 
larger values of q it is not possible to rigidize the lattice 
by randomly adding diagonals), and we suggest that the 
"fully-first-order" transition (i.e. a first-order backbone) 
only occurs in the special case of a perfect (undiluted) 
square lattice. 



IV. CONCLUSIONS 

We have compared three types of percolation transi- 
tion in two dimensions: the connectivity transition and; 
the generic rigidity transition on the triangular lattice 
and; the generic rigidity transition on the braced square 
lattice. A summary of our understanding is as follows: 
(i) The generic rigidity transition on triangular lattices 
is second order with i/ = 1.16 ± 0.03, < /? < 0.2, 
(3' = 0.25 ± 0.02 and; (u) The rigidity transition on 
the braced square net is first order with finite backbone, 
spanning cluster and cutting bond densities at the per- 
colation threshold. Only the value of /3 for the generic 
rigidity transition on triangular lattices remains contro- 
versial, due to the very strong finite size effects in that 
case. 

To illustrate the fact that our data is inconsistent with 
a first order backbone in the site-joint rigidity case, we 
have developed the following scaling argument. 

Assume that the backbone mass |Q scales as Mb '^ 
L^^ at Pc- If the backbone is compact then Db = d, the 
dimension of the system. The backbone mass is com- 
posed of red (or cutting) bonds plus "blobs" of over- 
constrained, or self-stressed bonds (See Fig. g). There- 
fore AIb — Mred + Mbiobs- Thc uumbcr of red bonds 
in thc backbone scales as Mred '^ L^^'', as analytical 
results [^ 40 1 and the simulations reported here show. 
Let us furthermore write Muobs = nuobs x rnuobs where 
nuobs is the number of blobs in the backbone, and muobs 
be the average number of bonds in a blob. Therefore 
L^^ ~ L^/^ + nuobs'muobs- Now, the AS backbone is an 
exactly isostatic body-bar structure, formed by rigid clus- 



ters (blobs) joined by bars (red bonds) so that count- 
ing of degrees of freedom is exact on it and so M^ed = 
iubiobs + 2ns — 3 I^Oj . Here ng is the number of sites in 
the backbone, that do not belong to a blob (see Fig. ||). 
This identity is known as Laman's condition [ p9| , p0[ , and 
results from the fact that each red bond acts as a bar 
and therefore restricts one degree of freedom, while each 
blob has three degrees of freedom and isolated sites have 
two. The backbone is a rigid cluster and therefore has 
three overall degrees of freedom. We do not need to 
know Hs for our argument. It is enough to notice that 
nuobs < Mred/S '^ L^^". We can thus write 



^Ds <L^/'^{l + mbiobs/3) 



(3) 



To this point, we have made no assumption about the 
character (compact or fractal) of the backbone, so it is 
valid in general. If the transition is second-order, there is 
a divergent length (for example, the size of rigid clusters), 
and we expect muobs to diverge with system size. There- 
fore a non-trivial value results for Db , as we find numer- 
ically. If on the other hand there is no diverging length 
in the system, then muobs — > constant for large systems 
and Db — d — 1/v exactly. We thus see that a compact 
backbone requires an extensive number of cutting bonds, 
and this in turn can only be satisfied if j^ = l/d exactly 
|42| . This is completely inconsistent with our data, and 
so the possibility of a first order backbone is remote in 
two dimensions. 

The first order rigidity transition exhibited by the 
braced square net seems to be atypical as we illustrated 
using a model which tunes continuously from that limit 
toward the generic triangular lattice. We found that even 
a small deviation from the braced square lattice limit 
leads to a behavior similar to that of the triangular lat- 
tice. It would be intriguing if there were a tricritical point 
at which first order rigidity ceases and second order rigid- 
ity sets in, but we have not found a model which exhibits 
that behavior. Nevertheless there are a large number of 
other rigidity models in two dimensions, so the possibility 
is not yet ruled out. 
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